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Abstract 

A collective Hamiltonian for the rotation-vibration motion of nuclei is considered, in 
which the axial quadrupole and octupole degrees of freedom are coupled through the cen- 
trifugal interaction. The potential of the system depends on the two deformation variables 
P2 and P3. The system is considered to oscillate between positive and negative /?3- values, by 
rounding an infinite potential core in the (/32, /33)-plane with (32 > 0. By assuming a coher- 
ent contribution of the quadrupole and octupole oscillation modes in the collective motion, 
the energy spectrum is derived in an explicit analytic form, providing specific parity shift 
effects. On this basis several possible ways in the evolution of quadrupole-octupole collec- 
tivity are outlined. A particular application of the model to the energy levels and electric 
transition probabilities in alternating parity spectra of the nuclei ^^°Nd, ^^^Sm, ^^"^Gd and 
i^^Dy is presented. 
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1 Introduction 



Shape deformations and surface oscillations in atomic nuclei determine from a geometric 
point of view the main features of nuclear collective dynamics jT]. The leading quadrupole 
mode manifests itself in all regions of collectivity providing vibrational, rotational, and 
transitional structures of the spectra. In addition, in some regions the manifestation of 
octupole degrees of freedom is superposed, leading to more complicated shape properties 
and parity effects in the spectrum of the system [2 • A variety of microscopic, geometric 
and algebraic model approaches have been applied in nuclear regions where the quadrupole 
and octupole degrees of freedom coexist [3]. 

In general, the problem of quadrupole-octupole collectivity is not easy to solve neither 
microscopically, mainly due to the breaking of reflection symmetry, nor geometrically, due to 
the difficulty in determining the total inertia tensor of the system. It is, however, simplified 
considerably if the axial symmetry is still preserved and if the octupole deformations are 
fixed appropriately with respect to the principal axes of the quadrupole shape. Further 
simplification is achieved if both degrees of freedom are separated adiabatically. It allows 
one to examine the manifestation of the octupole mode for fixed values of quadrupole 
parameters. In such a case the collective motion can be associated to the oscillations of 
the reflection asymmetric shape with respect to an octupole variable /?3 in a double-well 
potential [HE]- Then the parity shift effect observed in nuclear alternating parity bands 
can be explained as the result of the tunnelling through the potential barrier [HI [7j. This 
concept has been generalized for the case of simultaneously contributing quadrupole and 
octupole modes jH], as well as for the case of higher multipole degrees of freedom 0. In 
both cases the double-well potential was defined in terms of a variable carrying the relative 
contribution of the different degrees of freedom and not the absolute values of the respective 
deformation variables. In such a way the explicit form of the original potential in terms of 
the quadrupole P2 and octupole jS^ deformation variables was not given. As a consequence, 
some basic characteristics of the quadrupole and octupole modes and their interaction 
remain outside of consideration. Such is the behavior of the system in dependence on 
the quadrupole and octupole stiffness, as well as the limiting case of a frozen quadrupole 
variable. Another interesting question is, if and to what extent one may consider the 
presence of a tunnelling effect in the space of the octupole variable after the quadrupole 
coordinate (32 is let to vary. Some limiting cases in the shape evolution and the angular 
momentum properties of the system are also of interest in respect with the above. 

The purpose of the present work is to clarify the above questions by applying a simple 
explicit form of the collective energy potential as a function of the quadrupole and octupole 
axial deformation variables (52 and /^a. We examine the evolution of the potential shape in 
dependence on both degrees of freedom, as well as on the collective angular momentum. 
The geometric analysis suggests that the oscillations of the system in the two-dimensional 
case of simultaneous manifestation of the quadrupole and octupole modes are performed in 
a different way, compared to the one- dimensional case of a reflection asymmetric shape with 
a frozen quadrupole variable. We study the physical consequences of the two-dimensional 
oscillations and demonstrate their role in the rotation-vibration motion of the system. 

In particular, the explicit geometric analysis of the quadrupole-octupole potential sug- 



2 



gests a possibility for a coherent interplay between both collective modes. This allows the 
derivation of explicit analytic expressions for the energy levels and electromagnetic tran- 
sition probabilities applicable to nuclei in which an "equal" (coherent) manifestation of 
quadrupole and octupole degrees of freedom is considered. As a result, one is able to study 
in detail the respective effects in the structure of the spectrum. Below it will be shown 
that such a consideration can be applied reasonably to some nuclei in the rare earth region, 
such as the N = 90 isotones ^^°Nd, ^^^Sm, ^^^Gd and ^^^Dy. These nuclei are also a subject 
of interest fOJ from the point of view of the X(5) critical point symmetry ^2] between 
quadrupole vibrations [U(5)] and axial quadrupole deformation [SU(3)]. In the present work 
we shall, however, mainly consider the common quadrupole-octupole collective properties, 
which, in principle, can take place in various nuclear regions. 

In Sec. 2 the Hamiltonian of the coupled quadrupole and octupole modes is presented, 
together with the geometric analysis of the quadrupole-octupole potential. In Sec. 3 the 
Schrodinger equation is considered in the case of a coherent interplay between the two de- 
grees of freedom. The analytic solutions for several particular forms of the potential and 
the respective schematic spectra are given in Sec. 4. In addition, results of the model 
description of alternating parity spectra in ^^°Nd, ^^^Sm, ^^^Gd and ^^^Dy are presented. 
The electric transition probabilities are considered in Sec. 5, while in Sec. 6 a brief dis- 
cusssion of the influence of the 7 degree of freedom on the present results is given. Finally, 
a summary and concluding remarks are given in Sec. 7. 

2 Hamiltonian for the coupled quadrupole and oc- 
tupole modes 

We assume that the system is allowed to oscillate with respect to the quadrupole P2 and 
octupole axial deformation variables. In addition, both degrees of freedom are coupled 
through a centrifugal (rotation-vibration) interaction depending on the collective angular 
momentum /. The energy potential represents a two-dimensional surface determined by 
the variables P2 and P3. 

The quadrupole-octupole Hamiltonian describing the collective motion under the above 
assumptions has the form 

where the potential is 

;/(ft.A./) = ic.ft=4c3A= + 5^^||^. (2) 

with X{I) = I{I + l)/2. Here B2 and -B3 are the effective quadrupole and octupole mass 
parameters, and C2 and C3 are the stiffness parameters for the respective oscillation modes. 

The last term in Eq. Q provides a coupling between quadrupole and octupole degrees of 
freedom. Its denominator can be associated to the moment of inertia of an axially symmetric 
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Figure 1: (Color online) Schematic 3D plot of the potential U{P2, Ps, I), Eq- (0), in MeV for 
/ = 5 as a function of /?2 and P^. The parameter values are C2 = C3 = 100 MeV, ^2 = 0.1 
h'^MeV-\ and = 0.01 n'^MeV-\ 

quadrupole-octupole deformed shape, j/'Ce^arf+oci) _ ^j^^p^ _|_ QB^P^ 13 . Therefore, the 
constants ^2,(^3 > can be related to the mass parameters as ^2 = 3-02 and = 6B3. 
However, in the present study we do not impose this relation and below a more general 
correlation between d2, B2 and d^, -B3 is considered. The quantities d2 and d^ determine 
the contributions of the quadrupole and octupole modes, respectively, to the moment of 
inertia. Also, we remark that if the ground state of the system is considered (/ = 0), the 
potential 11(^2, P3, 1) should be taken by replacing X(J) do + X{I), with d^ being a 
constant. 

The Hamiltonian represents a two-dimensional generalization of the soft octupole 
oscillator Hamiltonian introduced in |14|, as well as of the one-dimensional octupole Hamil- 
tonian derived in [T3] . In the latter two approaches the quadrupole mode is assumed frozen 
as mentioned in Sec. 1. In this respect, Eq. corresponds to an extension in which 
the quadrupole coordinate is let to vary. Also, it corresponds to the quadrupole-octupole 
Hamiltonian in Refs. [HI E] • However, in the present work the potential energy (j2)) is taken 
in an explicit form depending on j32 and (3^ (including the harmonic oscillator part), while 
in Refs. jHl El a double-oscillator potential is defined in the space of polar coordinates. In 
this way, the explicit form of Eq. Q allows one to examine in detail the potential surface 
and its dependence on the model parameters and the collective angular momentum. 

Having in mind that the quadrupole deformation has the leading role in the rotation 
mode, we assume that its contribution to the moment of inertia is larger than the octupole 
contribution. This assumption corresponds to the condition d2 > d^, e.g. we can take 
d2 = 0.1 ^^MeV~^ and ds = 0.01 ^^MeV~^. Then, for comparable values of the deformation 
variables P2 and P3, the input of the quadrupole mode in the denominator of the centrifugal 
term will be larger than the octupole one. On the other hand, as it will be seen in Sec. 3, 
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Figure 2: Schematic numerical behavior of the potential U{j32, jSs, I), Eq. (j21), in MeV for 
J = 1 as a function of: (a) j32 at fixed /^a = 0.1; (b) at fixed /?2 = 0.25. The parameter 
values are = Cg = 100 MeV, d2 = 0.1 h'^MeV-\ and 4 = 0.01 h^MeV-\ 



this circumstance does not restrict the possibility of equal (coherent) contributions of both 
degrees of freedom in the mixed quadrupole-octupole oscillation mode. Moreover, since 
the rotation and vibration modes are coupled, the above condition might be not strictly 
imposed. In this meaning the considered values of d2 and c/3 provide only a schematic 
geometric analysis of the potential (j2I). 

Let us now examine the minimum of the potential energy in dependence on the model 
parameters. The set of extremum conditions for the coordinates of the two-dimensional 
minimum {p2min, Psmin) is 



_d_ 



U{P2,P3,I) 



and 



_d_ 



U{P2,P3,I) 



0; 



(3) 



d(3l 



UiP2,P3,I) 



> and 



UiP2,P3,I) 



> 0. 



It determines the following possible cases for the bottom of the potential, 

i) (33rmn = 0; /?2™n = ± [2X (/) / (rfsCa)] 

ii) /32™n = 0; P^rmn = ± [2X{I) / {d^C^)]" ' ' 

iii) /32mm 7^ and /^smm 7^ with the condition 

2X{I)d2 ^ 2X{I)d, 



and C?, 



^3miri) 



(4) 



(5) 



The shape of the potential corresponding to the case i) is illustrated in Fig. 1. It is 
characterized by two energy minima for > and /32 < separated by a well determined 
potential barrier. For given sign of (32 (we consider (32 > 0) the bottom of the potential is 
not separated in the P^- direction, allowing oscillations of the system between Ps > and 



Figure 3: (Color online) Schematic 3D plot of the potential 17(^2, Ps, I), Eq- in MeV 
for / = 5 as a function of (32 and jS^. The parameter values are C2 = 100 MeV, C3 = 10 
MeV d2 = 0.1 h'^MeV-\ and 4 = 0.01 h'^MeV-\ 



/^a < 0. This situation is illustrated in Fig. 2. We see that for a fixed physically typical 
value of (Fig. 2a) the barrier in the quadrupole space of j32 is very large. Thus it restricts 
the values of the quadrupole deformation within the half space P2 > 0. For a fixed typical 
(32- value (Fig.2b) the barrier in the octupole space of (3^ is relatively small. From Fig. 
1 it is seen that for some higher P2- values this barrier is reduced, and for P2 > P2min it 
disappears {Psmin = 0). 

In the case ii) the potential shape is the same as in Fig. 1, but the coordinates (32 and 
/?3 are exchanged. As far as the system is not considered to oscillate between positive and 
negative (32 deformations, this case is not of interest in the context of the present analysis. 

In the case iii) of non-zero P2min and Psmin, Eq. (0) imposes the relation 

TT = TT • (6) 

'-'2 <-^3 

It determines an elliptic form of the bottom of the two-dimensional potential surface given 
by (3l^J^2X{I)/{d2C2) + /5|™„/v^2X(/)/(rf3C3) = 1. The shape of the potential cor- 
responding to the case iii) is illustrated in Fig. 3. It suggests that the system moves in 
the two-dimensional space of the deformation variables (32 and (3^ by rounding the internal 
potential core. If a prolate quadrupole deformation is considered, the rounding is performed 
between positive and negative (3^ values in the space of (32 > 0. This situation can be 
considered as the two-dimensional extension of the one-dimensional case in which the (32 
coordinate is frozen. To explain this in detail, we consider a horizontal (equipotential) in- 
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Figure 4: (Color online) Schematic picture of the tunnelling and the rounding way in the 
/?3- variable. See Sec. 2 for further discussion. 



tersection of the shape in Fig. 3, which is illustrated schematically in Fig. 4. We see that if 
the quadrupole coordinate is fixed at some value of > 0, the motion in the octupole coor- 
dinate between positive and negative /^a-values is characterized by the tunnelling through a 
potential barrier (a vertical intersection of the core). When (32 is let to vary, the tunnelling 
is replaced by a motion along the curved way rounding the potential core. 

The above case iii) is of particular interest, due to the simultaneous presence of non- 
zero coordinates of the potential minimum in both degrees of freedom. It suggests that 
the oscillations in the quadrupole and octupole coordinates are involved in the collective 
motion on the same footing. As it will be seen below, such a situation appears to take place 
in certain nuclear regions. Moreover, it will be seen that the ellipsoidal symmetry in the 
potential bottom allows, under some additional conditions, a complete analytic determina- 
tion of the energy spectrum. This is why in the following we shall imply this case, unless 
something different is indicated. Also, we assume only the presence of prolate quadrupole 
deformations. This is why hereafter we consider only the /?2 > part of the space. 

Further, we examine the evolution of the potential shape with the angular momentum 
/. We consider the following two cases. 

I) The potential minimum (the 2-dimensional bottom) is allowed to change with / for 
fixed values of the stiffness parameters C2 and C3. 

II) The minimum is fixed, so that the values /32mm and (3^min determine an ellipse which 
does not change with the angular momentum. 

It is clear that in the case I) of fixed stiffness parameters, the quadrupole and the 
octupole deformations corresponding to the potential minimum should exhibit an overall 
increase in the denominator of Eqs. (0) with increasing I. 

In the case II) of fixed minima, the stiffness parameters C2 = C*2(-^) and C3 = 6*3(1) 
increase quadratically with I according to the right hand sides of (0). Then the substitution 
of Eqs. (0) into 0, leads to the following form of the quadrupole-octupole potential 

t/(/32,/33,/) = X{I) ((^^^/'+ /Ji^^^y + ^^^2 + ^3^2) • (7) 
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Figure 5: (Color online) 3D plot of the potential Ui{[32,l3z), Eq. (jHl), in MeV for / = 5 as 
a function of j32 and (3^ with (32min = 0.25, /^smm = 0.1, d2 = 0.1 /i^MeV"^, and da = 0.01 
^i^MeV-^ 



If the origin of the energy scale is fixed at the potential minimum, one has 

(8) 



Ui{P2,Pz) = f/(/52,/?3,/)-f/(/52™n,/93™n,/) 



{d2l3lmin + dAin)\d2l3l + ds/^f) 

We remark that the potential (|7j) includes the rotational contribution of the centrifugal 
term, which moves up the energy with increasing angular momentum /. On the other 
hand, in the potential (jHI) the explicit contribution of the rotational degree of freedom is 
diminished, so that the energy term keeps mainly the vibrational component. The shape 
of the potential ^ with /32 > is illustrated in Fig. 5. 

3 Model potentials and the Schrodinger equation in 
polar variables 

Further, it is convenient to introduce polar variables t] and cf) by taking 

P2 = I COS ; /?3 = sin , (9) 



'd2/d yds/d 
with d = {d2 + ds)/2. Considering P2 > 0, we have 




V = V ^2/^2 + dspi ; = arctan ( | , (10) 

where the "effective" deformation variable rj is defined with positive values 1] > 0, while the 
relative ("angular") variable is defined in the interval —tt/2 < (/) < tt/2. We remark that 
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the negative 0- values correspond to negative /?3. (The variable (3^ takes both positive and 
negative values.) 

Then the quadrupole-octupole Hamiltonian can be written in the form 
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Under the condition ©, the potential energy depends only on the effective deformation 
variable r] and on the angular momentum /, and not on the relative (angular) variable 0. 
Then in the case I) of fixed stiffness parameters one has 



(Case I), 



(12) 



where C is defined according to Eq. (0) as 1/C = ^2/(^6*2) = ds/dC^. 

In the case II) of fixed minima the potential term appears in the following two forms 
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X{I) 
X{I) 
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Immi 
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(Case II. A); 
(Case II.B), 



(13) 
(14) 



where Eq. (fT^ corresponds to the rotation dependent potential ((Tj), while Eq. (fT^ repre- 
sents the essentially vibrational term (jSI). The quantity T]min = (l/v^)('^2/52mm+'^3/^3mm)"^''^ 
is the value of the variable 77 in the potential minimum. In the following we shall refer to 
Eq. (fT^ as case II. A and to Eq. (fT^ as case II.B. 

Using the effective deformation variable rj we are also capable to examine a third case 
(III) of an infinite square well with an infinite core at zero defined as 
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(Case III), 



(15) 



where ?7w is a parameter determining the width of the well. 

Now we assume the following relation between the quadrupole and octupole mass and 
inertia parameters 
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(16) 



This leads to the following form of the model Hamiltonian 
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1^ 



(17) 



The assumption (fTB)) . which much simplifies the problem, suggests that d2 and are related 
to the mass parameters B2 and respectively, through the same coefficient d/B. By 
comparing Eq. (fTH|) and Eq. (jH)), we obtain C2/-B2 = C3/B3, or = ti;|, i.e. the assumption 
(fTB|) implies that both degrees of freedom, quadrupole and octupole, are characterized by 
equal angular frequencies uj2 and cus, respectively. This means that a coherent interplay 
between the two collective modes is assumed. In other words, the condition (fT^ suggests 
that the oscillations in the quadrupole and octupole coordinates are represented in the 
collective motion on the same footing. The quantity B in Eq. (jl6p has the meaning of the 
effective mass of the total quadrupole-octupole system. 

The Schrodinger equation for the Hamiltonian (fT7|) has the form 



2B r/2 



^2^2 Q 

drf ^ drj 



^{ri,^) + Ui{ri)^{ri, 



E^{ri, 



After dividing it by h'^/{2Brj^) and separating the variables i] and through $(77, 
■ip{ri)ip{(f)) we obtain the following two equations 



2B 



2B rj^ J 
52 

— <^(0) + k^ificj)) 



0; 
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(19) 
(20) 



where k is the separation quantum number. 



4 Analytic solutions and numerical results 

In the following we give analytic solutions of the above equations in the cases I-III with the 
potentials (IT^. E^ . (ITl|l . and (ITH|) . 

In the Case I, after introducing the potential ()12p into Eq. p9p we have 
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2B 
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25r/2 2 ' dr]^ 
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By introducing a reduced energy e = and a reduced angular momentum factor X(I) 



bX{I), with b 



2B 



we obtain Eq. (PT|) in the form 



1 d 



dr] 
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The effective potential appearing in the brackets of Eq. (|^^ is of a form similar to the 
Davidson potential [IS], which is analytically solvable [T71[T8J. Thus Eq. can be solved 
analytically and we obtain the following explicit expression for the energy spectrum 



En,k{I) = 



(23) 



where uj = yC/B and n = 0,1, 2,.... The eigenfunctions ip^rj) of Eq. (PT|) are obtained in 
terms of the Laguerre polynomials 



\ 



2r(n + 1) _a„2/2 s 2st2s/ 2\ 

— -, -e ''an L„ [ari ] 



(24) 



where a = VBC /h and s = ^k'^ + X{I)/2. 

Now we remark that Eq. (j2Up in the variable is solved under the periodic boundary 
condition (p{(j) + 27r) = (p{4>). On the other hand the assumption j32 > 0, which is equivalent 
to the consideration of an infinite potential wall at /32 = (or = ±tt/2), imposes the 
additional condition 



v{-n/2) = ^{n/2) = . (25) 

Eq. (PUI) has two different solutions satisfying the condition with positive, vr = (+), 
and negative, tt = (— ), parity as follows 

¥^+(0) = y27^cos(A;0) , A; = ±1,±3,±5,... ; (26) 
(^-(0) = .y27^sin(fc0) , = ±2,±4,±6,... . (27) 

Eq. (PU)) provides positive parity for the intrinsic wave function, while Eq. (fTfj) corre- 
sponds to a negative parity function. As a result, the intrinsic wave function appears in 
the form ^^{rj, 0) = ip{rj)ip'^{(f)). On the other hand, the TZV- symmetry of the total wave 
function of the system, ^> ~ <I>^(?7, 0)|/i^'M), has to be conserved. The TZ- symmetry of 
the rotation function \IKM) is characterized by the factor (—1)''^. For the total state of the 
system one has 7r(— 1)^ = 1. It follows that for J=even the quantum number k is allowed to 
take the values = 1, 3, 5, . . ., corresponding to the even function (j^ . while for J=odd 
one has = 2, 4, 6, . . . corresponding to the odd function ()27p. Thus, when the angular 
momentum is changed from J=odd to /=even and vice versa, the respective values of the 
quantum number k should switch between and This parity effect provides an 

odd-even staggering structure of the spectrum (j2Sl)- We consider that the lowest states of 
the system with respect to the variable are characterized by the lowest k- values, = 1 
for / =even and = 2 for / =odd. Therefore, the staggering behavior of the model 
spectrum is provided by the difference AA;^ = — /c^^^ = 3. 

In such a way the energy expression (j23p . with the parity-dependent quantum number 
fc, determines the structure of an alternating parity spectrum. The energy levels EQ^k{I), 
with n = 0, correspond to the yrast alternating parity sequence. The levels with n ^ 
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Figure 6: Schematic energy levels (a) and staggering pattern (b) for the spectrum ()23|) with 
n = 0. The parameter values are shown in part (a), oo is given in MeV/h, while b and do 
are in h'"^ and h'^ respectively. 



correspond to higher energy bands, in which the rotational states are built on quadrupole- 
octupole (mixed P2-P3) vibrations of the system. In this case, the states with even I appear 
similarly to the states of a higher /3- (quadrupole) band. Thus, the present model sug- 
gests that, in the nuclear regions with quadrupole-octupole collectivity, one may consider 
"octupole mixed" f3- band structures. We remark that, in the present model framework, 
the 7- bands are not included. This can be done in an extended formalism allowing the 
simultaneous consideration of the 7- variable. In addition, the octupole triaxiality can be 
taken into account. Then one may also discuss possible "octupole admixtures" in the 7- 
band structure. 

We can estimate analytically the staggering effect at higher angular momenta, where 
Ak^ <^ X{I). The square root term \Jk!^ + X{I) in Eq. can be expanded as \J X{I) + 

k'^ / {2\l X{I)). We see that the term k'^ / {2\J X{I)), which is responsible for the staggering 
effect, decreases nearly linearly with the angular momentum / since X{I) = bX{I) ~ 
1(1 + 1). We consider the quantity b = 2B/{h d) as a model parameter. The numerical 
behavior of the energy and the staggering effect for the spectrum is illustrated in Fig. 6. 
The staggering effect is illustrated in terms of the five point quantity 

Stg(/) = QAE{I) - AAE{I - 1) - AAE{I + 1) + AE{I + 2) + AE{I - 2) , (28) 

where AE{I) = E{I + 1) — E{I). The schematic staggering pattern suggests that the 
odd and even angular momentum sequences approach each other towards higher angular 
momenta. It outlines a trend for the forming of an octupole band. However, the linear 
decrement of the staggering amplitude is not enough to provide such a band structure at 
reasonable (experimentally observed) angular momenta. A similar situation is observed in 
rare earth nuclei, where the alternating parity levels approach each other without merging 
into a single band. 

On this basis, we applied Eq. ()23|) to describe the alternating parity spectra in the 
nuclei ^^°Nd, ^^^Sm, ^^^Gd, and ^^^Dy. The theoretical energies are obtained by taking 
En,k{I) = En,k{I) — -Sn,A;(0), with u = and X{I) = |[(io + /(/ + 1)], where the parameter 
do characterizes the potential shape in the ground state, as mentioned in the paragraph 
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Figure 7: (Color online) Theoretical and experimental energy levels for the alternating 
parity bands in ^^ONd (data from ^^^Sm (data from [20]), ^^^Gd (data from gT]) and 
^^^Dy (data from [22] )• The theoretical results are obtained by (PHj) with n = 0. The 
parameter units are as in Fig. 6. 



after Eq. The parameters lu, b and are adjusted to the energy levels by means of a 
least square minimization procedure. 

In Fig. 7 results for the energy levels of ^^''Nd, ^^^Sm, ^^^Gd and ^^^Dy are compared 
to the experimental data. The respective theoretical and experimental staggering patterns 
are compared in Fig. 8. In ^^°Nd [Fig. 7(a)], the levels with I = 9,11,13 are predicted. 
The respective staggering pattern for / > 5 [Fig. 8(a)] is also predicted. We see that in the 
nuclei ^^^Sm, ^^''Gd, and ^^^Dy the experimental patterns confirm the predicted behavior 
of alternating parity levels with increasing angular momentum. 

In the Case II. A, after introducing the potential fT^ into Eq. |T^. we have the equa- 
tion 



^^(,) + -^V^(,) + ^ 



E 



X{I) 



drf drj!, 



.Y(/),^^, 



min 



(29) 



which is solved in the same way as Eq. (PT|) of case I, and the respective energy levels are 
obtained in the form 



En,k{I) = 'h^ 



X{I) 



2 

vain 



2n + 1 + \/k^ + X{I) 



(30) 
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Figure 8: (Color online) Theoretical and experimental staggering patterns for the alternat- 
ing parity bands in ^^°Nd, ^^^Sm, ^^^Gd and ^^^Dy. The theoretical results are obtained by 
(j23|) with 77, = 0. The parameter units are as in Fig. 6. 



where n = 0, 1, 2, . . .. The eigenfunctions iP{t]) of Eq. ()29|) are of the same form as Eq. 



in case I, but with a = yX(J)/?7^-„. 

All considerations related to the 0- equation (pUj) and the quantum number k are the 
same as in case I. However, now we obtain a different behavior of the staggering am- 
plitude as a function of the angular momentum. This is seen after expanding the term 

\/ X{I)\J k"^ + X{I) of Eq. (jnni) in the form X{I) + k"^ /2. The appearance of the staggering 
effect is only due to the term k'^ /2. Since the difference A/c^ = 3 does not depend on J, the 
staggering effect will be characterized by a constant amplitude. The schematic behavior 
of the energy levels and the respective staggering pattern for the spectrum of Eq. (jHUj) are 
illustrated in Fig. 9. Indeed, we see from Fig. 9(b) that, after some slight increase in the 
beginning, towards the higher angular momenta the staggering amplitude saturates to a 
constant value. 

In Case II. B, the potential ()14p differs from the potential ()13|) of case II. A by the term 
— 2X(/)/((i?7^-^). The respective energy spectrum is 



2 

min 



2n + l + ^k^ + X{I)-^X{I) 



(31) 



with the wave function ip{ri) being the same as in case II. A. The expression ()31|) differs from 
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Figure 9: Schematic energy levels (a) and staggering pattern (b) for the spectrum of Eq. (jHUj) 
with n = 0. The parameter values are shown in part (a), B is given in h^/MeV, b in h''^, 
while rjmin is dimensionless. 




Eq. (j3(jp by the term —\JX{I) in the brackets. This term reduces the angular momentum 
dependence of the energy to a linear (vibrational) behavior. On the other hand, it does 
not affect the staggering effect. Therefore, similarly to the case II. A, the staggering pattern 
for the levels pip will be characterized by a constant amplitude. This is seen from the 
schematic numerical results illustrated in Fig. 10. 

In Case III [the square potential well (fTH|) ] the Schrodinger equation can be written in 
the form 



Id 

—ip{r]) + -—ipiv) + 



drj 



7] drj 



e — 



7]^ 



(32) 



where u"^ = k"^ + ^i^) and rj < rj^ . By introducing new variables through the definitions 



rjK and e 



we obtain Eq. (j32j) in the form of the Bessel equation 



Id 



dz 



z dz 



1 - 



ij{z) = 0. 



The spectrum of this equation is determined by the boundary condition tpuiVv, 
is given by 

.2 



K, 



(33) 
0, and 
(34) 
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Figure 11: Schematic energy levels (a) and staggering pattern (b) for the spectrum of 
Eq. ((23) with n = 0. The parameter values are shown in part (a), B is given in 
is in ^^MeV~^, while rj^ is dimensionless. 

where x^^n is the n-th zero of the Bessel function Jv{z). The eigenfunctions have the form 
4'u,n{v) = Cu,nJu{i^u,nV)j where Ci,^n are normahzation constants. The schematic behavior 
of the spectrum (jH^ and the respective staggering pattern are illustrated in Fig. 11. We 
remark that the staggering amplitude initially decreases, while towards higher / it saturates 
to a constant value. 



5 Electric transition probabilities 

The formalism developed so far allows the calculation of El, E2 and E3 transition prob- 
abilities for the energy spectra in the considered cases I-III. In the cases I and II, the 
reduced probability for an electric transition of multipolarity L from a state with angular 
momentum Jj to a state with If is given by 

B{EL-h^Ij) = -^ Y: \{^%,^M^vA)\M,{EL)\^lj^^X^,cl>))f , (35) 

where 



2T(n 



\\ r(n + 2s + l) 

The general form of the multipole operators Ai in the collective variables is given in 
The electric quadrupole and octupole transition operators for an axially symmetric 
nucleus are defined by the deformation variables /?2 and (3^ as 

M,{EL)=MlPlD^,{6), L = 2,3, (/i = -L, L), (37) 

while the El (dipole) transition operator is defined as pij-^T] 

M,{E1) = M,P,P,Dl^{e) , (/i = 0, ±1), (38) 
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where Mi {i = 1,2,3) are constants related to the respective intrinsic moments. In terms 
of the polar variables rj and (p the above transition operators read 

MAEl) = M/f^^£^Dl(e), (39) 

M,{E2) = M.^l^DliO), (40) 



d2/d 

M,m = M,^Dl{9) . (41) 
Ids/ a 



/'++ 


2 

7r . 


/ COS^ : 


(pd4> = 




2 

TT . 


r7r/2 

/ cos(/: 

l-n/2 


32 

.sin2(2</.)# = — , 


4r 


2 

TT . 


/ COS^ - 
I-tt/2 


16 

(/)sin0sin(20)(i0 = 

ISvr 




2 

TT . 


r7r/2 

/ cost/: 

l-n/2 


isin0sin(20)(i0 = -. 



In Eq. ()35|) the integration over the angles 6' involves an integral over three Wigner 
functions 2^, which leads to the Clebsch-Gordan coefficients (JjOLO|//0). The integration 
over the variable </> leads to the following constants 

(42) 
(43) 
(44) 
(45) 

The notations (++), ( ), and (H — ) correspond to the parities of the functions 

included in the integration. 

As a result, the reduced El, E2, and E3 transition probabilities between levels with 
|njJj) and \nflf) are given by the expressions 

B{E1,I, If) = h{WW\IfO)^S''{El,Ii ^ If), (46) 

B{EL, I, ^ If) = bL{WLO\IfOfS\EL, h ^ //) , (47) 
where _ 

POO 

S{El,k^If) = dv^if^{v)v'^i]{v), (48) 

POO 

S{EL,k^If)=i dr^i^lf^{r^)ri^,l,l^X^), (49) 







with L = 2,3. In Eq. (jl7|) the values of the integrals (021), (021), (03) are included in the 
constant h^. In Eq. (j46p the constant (j44j) is included in hi. We remark that if values of 
different kinds of transition probabilities are compared, or if branching ratios are considered, 
the quantities (jl^ - ljl^j) should be taken into account explicitly. 



17 




Figure 12: (Color online) Theoretical and experimental B(E2) transition probabilities (in 
W.u.) as functions of the angular momentum in the alternating parity spectra of ^^°Nd, 
^^^Gd and ^^^Dy. The data for ^^^Nd, ^^^Gd, and for Experiment [A] in ^^"^Dy are from [23; 
Experiment [B] in ^^^Dy is from [SH!- The theoretical results are obtained by Eq. (jTfj) . See 
Sec. 5 for further discussion. 

In the case of transitions between states of the yrast alternating parity band, |0/j) and 
\OIf) (with Hi = Uf = 0), we obtain the integrals ()48|1 and ()49|1 in the following simple 
analytic form 

QfT?^ T T\ 1 TjSj + Sf + 2) 

S{El,Ii If) = —— (50) 

«'^r(2s, + i)r(2sj + i) 

S{EL, J, If) = 4^^^=H^^ii4==. (51) 



r(2si + i)r(2s/ + i) 



where Si = il/2)^Jkf + X{Ii), Sf = {l/2)^k} + X(//), and a = VBC/h. 

In the case of the infinite square well potential (case III), the model wave function is of 
the form 
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Figure 13: (Color online) Theoretical and experimental B(E2) [part (a)] andB(El) [part (b)] 
transition probabilities (in W.u.) as functions of the angular momentum in the alternating 
parity spectrum of ^^^Sm (data from 123 )• The theoretical results are obtained by Eqs. 
and (gZI). 

while the integrals over the variable rj read 

roo 

S{El,Ii If) = drjJyXi^y,,n,vW'Juf{t^uf,nfV) (53) 

roo 

S{EL,Ii ^ //) = d7]J^X^y^,n,vWJuMuf,nfV)- (54) 

In general, the above formalism can be applied for a detailed analysis of the electric 
transition rates in spectra where the collective quadrupole octupole dynamics carries the 
characteristics outlined in the cases I-III of our study. In Figs. 12 and 13 we illustrate its 
application to E2 transition probabilities in the nuclei ^^°Nd, ^^^Sm, ^^^Gd, and ^^^Dy, as 
well as to the El transitions in ^^^Sm, in the framework of Case I. The results are obtained 
with the parameter sets given in Fig. 7. The quantity a = y/BC/h, appearing in Eqs. (fHUj) 
and (jHT|) . has been considered as a fitting parameter. The constant bi in Eq. ()4fij) has 
been determined so as to scale the theoretical El transition values with respect to the 
experimental data and takes the value bi = 1.2 x 10~^, while the constants 6^ in Eq. ^7\i 
have been set equal to 1. 

We see a good agreement between theory and experiment for the B(E2) values in ^^°Nd 
[Fig. 12(a)], i52Sm [Fig. 13(a)], and ^^^Gd [Fig. 12(b)]. In Fig. 12(c), the theoretical E2 
transition probabilities in ^^^Dy are compared to two different sets of experimental data, 
and [30^ (with no error bars reported in [301). that the theoretical values follow 

only the overall increase of the experimental data. We should, however, remark that the 
two sets of data diverge essentially, especially at the higher angular momenta. There is 
also some discrepancy between theory and experiment in the El transition values in ^^^Sm 
[Fig. 13(b)]. The results in Figs. 12(c) and 13(b) suggest that further examination of the 
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formalism, as well as of the experimental data, may be necessary. The further analysis of 
data on electric transitions in a wider range of nuclei will be the subject of future work. 



6 Influence of the 7 deformation mode on the (32-^3 
coUective motion 

As it has been mentioned in Sec. 4, the present model framework does not include the 7 
degree of freedom. Here we briefly discuss the possible ways in which this can be done and 
shortly estimate the influence of the 7 deformation mode on the collective motion in the 
space. The rotation energy of a system with presence of axial and triaxial quadrupole 
modes {P2 and 7 ) and axial octupole degree of freedom (P^) can be given by ^131 EI] 

1 ^ if 

Trot = 2 E (2) ' .(3) ' (55) 
1=1 -r di 

where 

Jf 452/5^in2 f 7 - ?V (^ = 1,2,3) (56) 



3 J 

are the mo ment-of- inertia components of the quadrupole shape about the axes 1, 2 and 3, 
while those of the axially symmetric octupole shape are 

jf ) = jf ) = QB.Pl ) = 0. (57) 

A simple estimation of the 7- influence can be done by assuming small variations of 
the system around 7 = 0, as in the case of the X(5) model [T^. Then the quadrupole 
moment-of-inertia components ()56p can be taken as 

jf^) = jm = 3B2(3l = 452/^2 sin^ 7. (58) 

As a result the rotation energy ()55|1 obtains the form 

I ( P - il /| \ 

= 2 \w,pi + iB,pi + arM^I ■ ^^^^ 

The first term in corresponds to the centrifugal term in the quadrupole-octupole po- 
tential (0). The second term in provides the influence of the 7 mode on the potential. 
After taking into account Eq. ()59|) with ^2 = 3i?2 and d-^ = 6-B3, and by including a 7- 
oscillation term, the potential can be generalized in the form 



2 2 2 d2P2 + dsp^ 4ct2P2 sm 7 



where K is the projection of the angular momentum on the body-fixed z-axis. Then for a 
fixed value of 7 the extremum conditions (3) and (4) in Sec. 2 provide the following cases 
for the bottom of the potential (lUU)) in the P2-P3 space. 
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1) Psmin = with 



2min 



± 



d2C2 



^ ^ ' 2 V4sin2 7 



1/4 



- 1 



(61) 



2) /?2min 7^ and /^smin 7^ with the condition 

[2X(/)-ir^K I 3ir^ and r - " -^^]^3 ..^l 

The following comments can be done on the above result. 

i) The appearance of /3| in the denominator of the second term in ()60|) divides the (32-^3 
space into two half-spaces, /?2 > and P2 < 0, separated by an infinite potential barrier 
at P2 = 0. For this reason the potential minimum with /?2min = and /?3min 7^ does not 
appear. 

ii) Eqs. ()61|) and (jU^ illustrate the ways in which the term involving the 7 deformation 
mode can shift the position of the potential minima in the /52-/33 space (compare with cases 
i)-iii) including Eq. (jSj) in Sec. 2). Note that for K = the influence of the 7 mode on 
the P2-P3 potential shape automatically disappears. This is a limit in which the P2 and 7 
degrees of freedom are weakly coupled and can be adiabatically separated, which is implied 
in the framework of the present work. The involvement of the K = 2 configurations in the 
collective motion implies the consideration of a strong /?2-7 coupling |i32] . 

iii) The involvement of the 7 degree of freedom in the above way would influence the 
correlation between the axial /?2 and (3^ variables due to the appearance of the second term 
in C2 [Eq. (jU^ ]. as a consequence of the last term in the potential (jUUj) . Now, in terms of 
the polar variables, the potential will depend on both rj and (j), so that the variables in the 
Schrodinger equation cannot be directly separated. This could be done in a way similar to 
the adiabatic separation of the f5 and 7 degrees of freedom in the X(5) model framework 
|12j . as well as in the framework of the AQOA model [S3]- Alternatively, the problem could 
be solved numerically in a way similar to the approach of Ref . [22] • 

A more general way to examine the influence of the 7 deformation mode on the qua- 
drupole-octupole motion of the system could be based on the complete form ()56|) of the 
quadrupole moment-of-inertia components, so that the 7 variable would not be limited in 
the vicinity of zero. Furthermore, non-axiality of the octupole degree of freedom can be 
considered. Any efforts in these directions should be based on numerical solution of the 
problem. 



7 Summary and conclusions 

The present study outlines some dynamical properties of a system with simultaneously 
manifesting quadrupole and octupole degrees of freedom. We remark that the obtained 
results represent a restricted class of exact analytic solutions of the problem. This is due 
to the correlation (jl6p between the mass and the inertial parameters, which essentially 
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simplifies the Hamiltonian (fTT|) in the form of Eq. (fT7|) . In addition, the correlation (jH)) 
between the inertial and oscillator parameters brings the potential in a form depending on 
the "effective deformation" variable rj only, and not on the relative "angular" variable 0, 
thus allowing an exact separation of variables in the Schrodinger equation. As it is explained 
in Sec. 4, the above correlations provide a coherent interplay between the quadrupole and 
octupole collective modes. In this respect, the presently considered potentials suggest some 
specific properties of quadrupole-octupole collectivity, which can be developed in various 
nuclear regions. 

However, despite the above limitations (the necessary price we pay for solving the prob- 
lem exactly) we were able to identify a region of nuclei where the assumed "equal" presence 
of quadrupole and octupole degrees of freedom can take place in the collective motion. We 
found that the structure of the spectrum in the case of the potential (fT^ . illustrated in 
Fig. 6, is similar to the structure of alternating parity bands in some rare earth nuclei. 
On this basis, we have reproduced quite accurately the energy levels and the staggering 
patterns in the nuclei ^^°Nd, ^^^Sm, ^^"^Gd, and ^^^Dy, as demonstrated in Figs. 7 and 8. 
In these spectra the reduction of the staggering amplitude indicates the trend of forming 
octupole deformations towards the higher angular momenta. However, the slow decrease of 
the parity effect does not allow this to happen at reasonable (observed) angular momenta. 
The B(E2) transition probabilities have been also described with a reasonable accuracy 
(Figs. 12 and 13), while the result for B(E1) transition probabilities in ^^^Sm [Fig. 13(b)] 
suggests further tests of the formalism and analysis of additional experimental data. 

We remark that the energy expression ()23p cannot reproduce the complicated beat 
staggering effects observed in the octupole bands of light actinide nuclei . The latter have 
been described [H], with good accuracy, by the use of the Quadrupole-Octupole Rotation 
Model Thus, compared to the light actinide region, the application of expression 
indicates a different behavior of the quadrupole-octupole collectivity in the rare earth nuclei, 
with a less developed octupole deformation and a more strongly pronounced octupole vibration 
mode. The present analysis suggests that in this case the coherent (equal) contribution of 
quadrupole and octupole oscillations can take place in the collective motion of nuclei. 

The potentials with fixed energy minima (cases II and HI) can be related generally to 
a situation in which the vibrational and rotational degrees of freedom are weakly coupled. 
Then the rotational angular momentum slightly affects the quadrupole-octupole vibration 
motion, which suggests a constant (or nearly constant) behavior of the staggering ampli- 
tude. In particular this is well seen in the case II. B, where the direct contribution of the 
rotational motion is excluded. Thus, the spectrum illustrated in Fig. 10 suggests an essen- 
tially quadrupole-octupole vibrational motion of the system. Cases II. A and HI, in which 
the rotational mode is taken into account, suggest quadrupole-octupole vibrations with 
an adiabatically manifested rotational motion. We remark that the constant staggering 
patterns, illustrated in Figs. 9 and 10, are in some meaning idealistic C3iSGS, clS far as the 
current experimental data do not show such a strong persistence of the parity effect at high 
angular momenta. On the other hand, the square well potential of case III appears to be 
applicable to examining the possible critical behavior of the quadrupole-octupole collectiv- 
ity in different nuclear regions. Studies in this direction have been implemented recently 
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in the light actinide nuclei [221 • We suggest that further analysis of experimental data for 
quadrupole-octupole spectra would be of use for testing the prediction of the staggering 
pattern illustrated in Fig. 11 for the case of the square well potential (fT3j) . 

It is important to note that the present exactly solvable model can be naturally ex- 
tended, beyond the "coherent interplay" assumption, to a more general non-analytic prob- 
lem in the following two ways. First, we can release the correlation (jH)) between the inertial 
and oscillator parameters, allowing the potential to depend on the variable 0. Then the 
problem can still be transformed into a form having an analytical solution, by performing 
an "approximate" separation of variables, as done in and in the framework of the X(5) 
symmetry model [Ej- The second extension would be to release the correlation (fTI)|) be- 
tween the mass and inertial parameters. This would allow us to examine different ways in 
which the coupled quadrupole and octupole degrees of freedom enter the collective motion. 
In this case, however, more sophisticated mathematical and numerical techniques have to 
be sought to solve the problem. 

Finally, we also remark that the developed formalism contains several limits. Thus, when 
the quadrupole variable is frozen to some stable quadrupole deformation, the potentials 
and the spectra of the cases I-III transform to the respective ones appearing in the one- 
dimensional problem ^j. Another interesting limit can be obtained by appropriate 
parameter values, for which the difference Ak"^ is negligible compared to X{I) for all angular 
momentum values. Then the staggering effect vanishes, and the odd and even angular 
momentum sequences appear in a single non-perturbed collective band. For example, if 
such a transition is performed in case III, the spectrum presented in Fig. 11 is reduced 
to the structure supposed to correspond to the transition between octupole vibrations and 
stable octupole deformation, in which a single octupole band is formed [2Sj. It is also of 
interest to take into account the non-axiality of the quadrupole and/or the octupole degree 
of freedom, in order to examine how the present results are modified. Studies in these 
directions will be the subject of further work. 
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